Propagation of sound in two species Bose-Einstein condensates 
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We investigate the propagation of zero sound in two-species interpenetrating Bose-Einstein con- 
densates. In very elongated clouds, this propagation is shown to be essentially one-dimensional. We 
also present ID numerical experiments that exhibit the existence of two different sound pulses with 
velocities consistent with the theoretical predictions. 



The dynamical behaviour of one-component Bose-Einstein condensates (BEC) has been the subject of numerous 
experimental and theoretical studies In particular, measurements of the velocity of sound pulses propagating in 

a cloud of 23 N a atoms have been shown to agree with the prediction of Bogoliubov theory 0. On the other hand, 
Myatt et al. have reported the observation of a double condensate created by thermal contact between two different 
hyperfine states of 87 Rb [Q. This technique of sympathetic cooling combined with the design of an optical dipolc 
trap H opens up the way to the realization of different multi-component condensates. Several theoretical works have 
already been devoted to the study of the ground state properties of two-species Bose-Einstein condensates (TBEC) 
H as well as to the determination of the frequency of the lowest collective excitations which correspond to global 
oscillations of the clouds j7|-|| . 

In this paper, we study the propagation of zero sound in large two species condensates. Such systems exhibit at 
long wavelengths two phonon branches. As in the case of pure condensates [ ^0|JTl| ], we show that in long cigar-shaped 
clouds, this propagation takes place essentially along the long axis. The sound velocities then take a simple form that 
lends itself easily to numerical and experimental verifications. 

Ultracold and dilute TBEC are well described by coupled Gross-Pitaevskii (CP) equations for the macroscopic 
wavefunctions: 
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together with the normalization conditions 



dr\^ t \ 2 =N u (3) 

where the numbers of particles are conserved. 

In these systems, the low energy binary collisions can be characterized by the s-wave scattering lengths ay. The 
interactions strengths are given by 

Uij =2-Kh 2 aij [— + — ], (4) 
rrii rrij 

Uij is positive when the interaction is repulsive and negative for attractive potentials. Few works have been devoted 
to the determination of the scattering lengths between unlike alcalis. In the following, we restrict ourselves to the 
case Uu > U22 > 0. The second terms in Eqs. (Q) and (||) are the anisotropic harmonic confining potentials: 
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where ^^"j- 1 = ^ \ 12 1 ; gj is the Lande 3 factor. We only consider mixtures of alkali atoms for which g\ = (72 and Eq. 
(5) then becomes 

V e (r) = ?-R 2 ; R 2 = x 2 +y 2 + X 2 z 2 , (6) 



where k = rriiU 2 ^. 

As in the MIT experiments ]^), we investigate the case of very anisotropic cigar-shaped potentials for which the 
asymmetry parameter A = is very small. The long axis lies in the z direction. 

In order to discuss the behaviour of collective excitations, it is convenient to introduce the Madelung transformation: 

Vi^yfce***, (7) 

where n, is the number density of the species i and, in absence of vortices, the phases cf>i play the role of potential for 
the two irrotational velocities v t — -^-V_4>i- In terms of these variables, the coupled GP equations take the following 
simple form: 

^+Y.(n^) = 0, (8) 



nVe(D + ]T U lf i 3 + ^ + p qi } = 0, (9) 
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Pqi = - — ^. V 2 (V^) (10) 



where 



is the quantum pressure. 

When the numbers of atoms are sufficiently large, the density profiles become smooth and the kinetic energy is 
small with respect to the interaction and potential energies. It only takes a significant contribution at the boundaries. 
When this term is neglected (p q i — 0), one can obtain simple analytical expressions for the density distributions in 
the atomic clouds (Thomas- Fermi approximation); in the case where U\i > U22 we get 

n T 2(r)='^[R 2 -R 2 }; forR<R 2 , (11) 



and similarly for n\ 



and 



n T ,{r) = ^[Ri~Rl] + ^[Rl-R\ forR<R 2 , (12) 



riTiir) = — — [Rl-R 2 ]; farR 2 <R<Ri, (13) 

and the profiles vanish elsewhere. 

In these equations, 5Ua = Ua — U12 and A = U\\Uii — U 2 2 . The radii R\ and R 2 of the clouds in the transverse 
plane at z = can be determined with the use of Eq. (||) . 

In free space (V e = 0), the binary mixture is thermodynamically stable with respect to phase separation if 

UuUaa ~ U 2 2 > 0; A > 0. (14) 

In the presence of the confining potential the two coupled clouds interpenetrate like miscible fluids as long as 
5U 22 > and the species with the largest self-interaction term (eg. 1) is more spread (cf. the bottom image in Fig. 
(1)). FromEq. (|l|) one can see that the species 1 experience an inverted harmonic potential near the origin (i? < R 2 ) 
when SU 22 < 0. The peak density of 1 does not coincide anymore with that of the other species, but it takes place at 



2 



R = R 2 . If the self- interaction 5U 22 is further decreased the species 1 is expelled from the center of the trap and form 
a shell at the periphery jsUli}]. In the following, we consider the case of two interpenetrating condensates, i.e. A > 
and 5U 22 > 0. 

It is useful to first recall the main formulas describing the propagation of sound in absence of external potential. 
The densities of the uniform condensate are then equal to the peak densities at r — obtained from Eqs. jTl| ) 
and jl^). The frequencies lu k of the perturbations about this uniform condensate and proportional to g^^-n-rl are 
given by the following dispersion relations: 



u; K± = -{[Un + U 22 ] ± [(Un - U 22 ) 2 + W^] 1 ' 2 } 1 ' 2 , (15) 
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Taking U X2 = and U n = or U 22 = in Eq. @, one recovers the standard Bogoliubov dispersion relations 
for the pure condensates of each species. In a TBEC, the mode corresponding to composition inhomogeneities is 



propagative and not diffusive as in the case of normal binary mixtures. Both frequencies in Eq. (15) are gapless and 
in the long wavelength limit (k — * 0), they yield two phonon-like branches: 

u>±(k — > 0) — > c±k, (16) 

where the corresponding sound velocities are given by Jig] 

c± = —7={{coi +c 2 ) ± (c i - c 2 ) H J ' I , (17) 

V2 mim2 

where coi = ( Ui ^ m ) 1 ^ 2 is the Bogoliubov expression for the sound velocity of a pure condensate of the species i. 
These collisionless modes correspond to the Goldstone modes which proceed from the spontaneous breakdown of 
the gauge symmetries — > e***$i in the composite condensate. In the mode associated to the frequency uj + , the 
two components interfere constructively whereas this interference is destructive in the case of As a result, the 
corresponding sound velocity c_ can become imaginary signalling a long wavelength instability leading to phase 



separation when A < |lq |. This condition is consistent with the thermodynamic criterion (cf. Eq. (14)). 

In the presence of the confining potential, the linearized equation for the perturbations <5rii(r, t) about the stationary 
profiles n e i take the following form after elimination of the velocity fields: 

d 2 5n 1 Marie! U 12 n el 

— -5- = V.{ Vdni} + V.{ Vdn 2 }, (18) 

at mi mi 



d 2 5n 2 Mi 2 n e2 M 22 n e2 

a-t9 = V.{ Vdni} + V.{ Vpn 2 \. (19 

ot 1 m 2 m 2 

When the densities are sufficiently high, the profiles n e i are well described by the Thomas-Fermi approximation (Eqs. 
(|IT|)-(|l3|)). This dynamics is greatly simplified in very elongated traps such as those used in the MIT experiments on 
pure BECs ||. Since the transverse dimensions are smaller than the longitudinal ones, the density profiles across the 
trap quickly reach their equilibrium values. After this short transient, the propagation is essentially one-dimensional 
and it can be described by a local density perturbation that depends only on the longitudinal dimension 5rii(z;t) 
ijniyrill . Moreover, at the lowest order in the smallness parameter A, one may neglect the variation of the equilibrium 
profiles n ei in the z direction (cylindrical geomety). Eqs. (jll]) and ( |l9| ) then become after integration over the cross 
section of the cloud of species 1 perpendicular to the direction of propagation: 

d 2 5ni _ Un <n e i > d 2 Srii Ui 2 < n e i > d 2 Sn 2 
Ot 2 ~ mi [ 9z 2 j mi 1 dz 2 h 1 ' 
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d 2 6n 2 _ Ui2<n e2 > d 2 6n 1 U 22 < n e2 > d 2 Sn 2 . , . 

dt 2 ~ m 2 [ dz 2 )+ m 2 [ dz 2 h 1 ' 

where 

2 f Rl 

Uei = ~n2 \ rdrneiir). (22) 
H i Jo 

As a result, sound propagation in such very anisotropic systems can be described by coupled one-dimensional wave 
equations. The corresponding sound velocities are still given by Eq. ([L7|), but the densities of the uniform condensates 
are now replaced by the average densities as over the large transverse section of the cloud. 

The effect of a density inhomogeneity on such two-species Bose-Einstein condensates can thus be captured by 
one-dimensional numerical experiments. The initial conditions were created by imaginary time integration of the 
coupled GP equations in the presence of an harmonic confining potential. Once the initial condition was obtained, 
we perturbated the condensate thus causing the appearance of inhomogeneities in the center of our mixtures. This 
perturbation was created by instantaneously adding to the initial external potential a Gaussian tip in the middle of 
the confining region. It mimics the laser beam that is focused in the experiments in the center of the trap to create 
a repulsive dipole force Q . A variant of this experiment consists of creating the initial condensate with an external 
potential that includes the Gaussian barrier and eliminating it instantaneously at time t = 0. Both constructions give 
qualitatively the same results. We show one such integration on Fig. (1). The parameters were chosen to guarantee 
that the mixture is thermodynamically stable and the clouds interpenetrate. For these parameters, the Thomas-Fermi 
approximation permits to calculate the extension of the condensates which are 462 and 570 fim. The numerically 
created condensates have widths of approximately 450 and 550 /xm, indicating a fairly good agreement between both 
aproaches. At t — we have switched on the above mentioned Gaussian barrier to perturbate the system. This 
generates two pulses which propagate symmetrically outwards. In the beginning one sees only one hump (on each side 
of the mixture). However, as time increases, one can distinguish two humps one of which outdistances the other. The 
position of these density minima varies linearly with time and the speeds of propagation can thus be extracted from 
the simulations. They are of approximately 10.4 and 8.1 cm/s. When we calculate them by using the Thomas-Fermi 
peak densities into the expressions for the sound speed, we find 11.7 and 9.4 cm/s. Again there is a fairly good 
agreement between the integration of the coupled GP equations and Thomas-Fermi approximation results. 
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Figure 1. Numerical integration of the coupled GP equations. The bottom image corresponds to the initial condition. 
The upper images show the development in time (increasing upwards and given in /xs) after the perturbation has 
been applied. On the y-axis we plot the density of these condensates in ,uto~ 3 , while the abscissa represents their 
spatial extension in fim. The dashed line corresponds to condensate of particles 1. The dot-dashed line corresponds 
to condensate of particles 2. The sum of the densities is plotted as a solid line. For this figure the masses of both 
species are 10~ 2e Kg. Each species contributes with 5 10 6 particles, and the frequency of the harmonic oscillator is 
equal to ZQQHz. Particles of type 1 and 2 present s-wave scattering length of the 1.5 10~ 3 and 8.0 10~ 4 /itm, while 
the interspecies length is of 1.0 10~ 4 /im. At t = a Gaussian barrier was switched on to perturbate the mixture. 
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